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We investigate the Affleck-Dine mechanism when multiple flat directions have large values simul- 
taneously. We consider in detail the case when both LHu and HuHd flat directions are operative 
with a non-renormalizable superpotential. In case Hubble induced A-terms are present for these 
two flat directions, their initial values are determined completely by the potential and there are no 
QQ , ambiguities how they are mixed. Moreover, CP is violated even when the Hubble parameter is large 

• due to the Hubble induced A-term and cross coupling in F-term, so that the lepton asymmetry is 

' generated just after the end of inflation. As a result, compared with the case of single flat direction, 

^SJ , the resultant lepton-to-entropy ratio is enhanced by a factor of Hosclm-^j^, where i/osc is the Hubble 

, ^ . parameter at the onset of oscillation and m3/2 is the gravitino mass. However, when Hubble induced 

' A-terms do not exist, there remains indeflniteness of initial phases and CP is violated spontaneously 

^ by the phase difference between initial phase and potential minima of the hidden-sector induced 

. A-terms. Therefore, CP-violation is not effective until the onset of the oscillation of scalar flelds 

y—^ • around the origin and there is suppression factor from thermal effect as is the case of single flat 

I direction. In this case, the amplitude of baryon isocurvature perturbation imposes constraints on 

, the model parameters. 

i-^ ' PACS numbers: 

53^' I. INTRODUCTION 

I The origin of the observed baryon asymmetry of the Universe [l| is one of the most important problems in both 
^— I . cosmology and particle physics which cannot be explained within the standard model of particle physics (SM) 0]. 
^ ' Inflation in the early universe , which is required in order to solve many cosmological problems such as the flatness 
and the horizon problems and to account for the origin of primordial fluctuations, washes out any preexist baryon 
■ asymmetry, so baryogenesis [1, Q must take place after inflation but before the big bang nucleosynthesis (BBN) Q 
epoch. 

As for the physics beyond SM, one of the best motivated models is supersymmetry (SUSY) Jj], because it can solve 
] the hierarchy problem and realizes the gauge coupling uniflcation. Generic supersymmetric models have flat directions 
> in fleld space where scalar potential vanishes at the classical level and their dynamics is very important in the broad 
' context. The Afflcck-Dinc (AD) mechanism @, H, [l^l of baryo(lepto)genesis makes use of scalar flelds having a large 
• • , amplitude along the flat directions. In this scenario, the angular momentum, or the rotation around the origin, of scalar 
. ^ ' fields that carry baryon or lepton number represents the number density of the field. Note that lepton asymmetry 
^> , can be converted to baryon asymmetry through the sphaleron effect [ll[, which violates B + L sd the electroweak 
5^ ' scale, where B and L are baryon and lepton charges, respectively. Therefore the source of baryon asymmetry in this 
scenario is large values of the scalar fields as the initial conditions and the phase difference between the potential 
minimum and the initial condition, which represents CP-violation. Later Dine et al.p^ analyzed behaviors of the 
field that parameterizes a flat direction (We call it the AD field.) during and after inflation taking non-renormalizable 
superpotential of the form 

n+3 

and SUSY breaking in the early universe into account. Here (p is the AD field and M is some cut-off scale. In this 
scenario, SUSY breaking in the early universe induces both a possibly negative mass term of the form, —cH'^\ip\'^ , 
and CP-violating A-terms of the form, {aH / M)Lp'^^^ + h.c. where c (real) and a (complex) are numerical factors of 
order unity and H is the Hubble parameter. The negative mass term assures the flat direction to have a nonvanishing 
vacuum expectation value (vev) and the difference of phase between a and the coupling constant of A-term from 
hidden sector of the form, {Amz/2 1 M)ip'^^^ + h.c, which is a consequence of SUSY breaking at low energy, provides 
CP-violation. When the Hubble induced A-terms do not exist, initial value of AD fleld in the angular direction is 
determined randomly during inflation and CP is violated spontaneously by the phase difference between initial phase 
of AD fleld and phase minima of the A-term from hidden sector. In this scenario, in both cases, the net baryon-to- 
entropy ratio is in proportion to the cut-off scale and the reheating temperature. Furthermore, it has been pointed out 
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that thermal effect induces an early coherent oscillation of AD field and suppresses the resultant baryon asymmetry 

[3, [HQ [a. 

Those analyses restricted only in the case where the configuration of flat direction can be parameterized in terms 
of one complex scalar field. There are, however, many flat directions even in the minimal supersymmetric standard 
model (MSSM) [H], some of which carry B — L charge but others do not. Therefore it is very important which flat 
direction, if any, is selected as the AD field. Moreover, there are multiple flat directions which do not give rise to any 
F-term in the renormalizable limit ^16*1. Such directions can get large values at the same time and we can no longer 
parameterize them by one scalar field. If some of them carry B — L charge and others do not, the degree of the mixing 
of multiple flat directions directly affects the net baryon asymmetry. It is not trivial whether the simple one-field 
analysis is applicable in such a case. 

The case where m ultip le flat directions have large values in the MSSM and its extensions have been considered 
in the literature [13, llSl [l9l . |20| . In this paper, we improve the evaluation of the potential and clarify how final 
lepton asymmetry is determined in the multifield context. In particular, we highlight the difference in the origin 
of CP-violation between the single-field case and the multiple-field case. We also study how much pre-inflationary 
information survives after inflation. 

As mentioned above the general structure of the flat directions of the scalar potential in MSSM has been classified 
by Gherghetta, Kolda and Martin 16], but not all these flat directions remain flat in cosmological setting where 
scalar fields are expected to be distributed rather randomly [21}. For example, if Higgs fields have large values, only 
LiH^[i — 1,2,3) and H^H^i flat directions remain flat and the other fields are driven to vanishingly small values. 
Here Li is left-handed slepton multiplet (z is family index), is the Higgs multiplet that couples to u-type quarks 
and Hfi is the Higgs multiplet that couples to d-type quarks and leptons. In this paper, we concentrate on such a 
situation, which is similar to the case [Tj]- We investigate the potential for these scalar fields in detail and estimate 
the net lepton asymmetry. We find that Hubble induced A-terms play critical roles in both the survival of pre- 
infiationary information and the dynamics of AD fields. If Hubble induced A-terms exist, our result for the resultant 
lepton asymmetry is consistent with the case of [13], but the dynamics of one of the scalar fields is slightly different. 
Moreover, if Hubble induced A-terms do not exist, there are almost no differences to the dynamics of scalar fields 
between single flat direction case and multiple flat directions case. In such a case, we must pay attention to the baryon 
isocurvature perturbation as is the case with Ref. [23] 

The structure of the paper is as follows. In the following section we present the relevant part of the MSSM in which 
we are interested. Then it is argued how CP-violating term is introduced and what is responsible for CP-violation. 
Furthermore we show how multiple flat directions are mixed in our model and how natural we consider the case 
where multiple fiat directions have large value. In section Hill AD mechanism in our model is discussed. The classical 
evolution along two fiat directions is studied. In section IIVI we comment on the case when the Hubble induced 
A-terms do not exist. We also discuss the baryon isocurvature perturbation. Finally, section |V] is devoted to the 
conclusion of the dynamics of scalar fields and leptogenesis in our model. 



II. MODEL 
A. Potential 

We adopt a non-renormalizable superpotential of the form, 

= ^{UH^){L,H^) + hE.[HM{HuHd), (2) 

in addition to the lepton sector of the renormalizable superpotential of MSSM [3], 

W^cn^eiYeijLjHd, (3) 

where e is right-handed slepton and ye is the Yukawa-coupling for lepton sector. Here M is some cut-off scale, \Lij 
and \h are coupling constants and i is family index. For example, if there is a right-handed Majorana neutrino 
multiplet whose mass M is heavier than Hubble parameter, we can acquire the first term of ^ by integrating out 
the right-handed sneutrino. We choose a basis which Xlij is diagonal. Note that this superpotential gives a neutrino 
mass when one of the Higgs fields gets a vev. 



(4) 
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Here = + (ff0)2 ~ 174GeV and tan/3 = {H^)/{H^). (0 means the neutral component of Higgs fields.) To 

realize inflation, we introduce an inflaton sector besides the MSSM sector. Couplings between the inflaton sector and 
the MSSM sector arises only with the gravitational strength in the supergravity scalar potential, 



V 



and assuming non-minimal Kahler potential, 



5Ki 



and ^i^2 - 7^/|0a|'+h.c., 



(5) 

(6) 
(7) 



in addition to the canonical terms. Here (f)a is scalar fields of MSSM and / is the inflaton, a is a label to distinguish 
the scalar fields and we assume the F-term of inflaton dominates the energy density during inflation. Note that while 
5Ki cannot be forbidden by any symmetries, some symmetries of the inflaton forbid 5K2. The resultant potential for 
Li , Hu and during inflationary era is 



IT 



UK 



^HanHjHl +h.c. 



(8) 



where Ca's are real and and an are complex parameters with their magnitude presumably of order of unity. Terms 
including H arise due to the coupling of flat direction to the F-term of the inflaton Fj. If 6X2 is forbidden, the Hubble 
induced A-terms, 



2M 



HaLLjHl + h.c. 



Ah 
2M 



HaHHjK + h.c. 



(9) 



vanish. Hereafter, therefore, we consider the two cases with and without Hubble induced A-terms. In the remainder 
of this section and section Hill we examine the case when the Hubble induced A-terms exist. In section HVl we discuss 
the influence of existence or nonexistence of Hubble induced A-terms to the dynamics of the scalar fields and the 
resultant lepton asymmetry. 

Soft terms from hidden sector SUSY breaking, 



^ms,2AL^Hl + h.i 



^m,/2AHKm + h.c. 



(10) 



exist at the same time. Here rria and 1713/2 are masses of scalar fields and gravitino respectively, while Al and Ah are 
complex numerical factors with their magnitude of order of unity. Scalar potential also includes D-terms, which are 



Vd 



(D^Y+Dl+Dj+Dj), 



Di^^l E(^r^i- + + (Hi* Hi + Hi* Hi) + {Hi* Hi + i/]*i/]) 

i J 

Y,iLl*L\ - Li*Ll) + (Hi* Hi - Hi* Hi) + [h]* Hi Hi* hI) 

i J 

{ E(|Lll^ - \L\\') + i\Hif ~ \Hif) + {\Hif iHif)] , 



D 



9i 
2 

9^ 
2 



(11) 
(12) 

(13) 

(14) 

(15) 
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in this case. Here g and g' are the gauge couphng constants and superscripts t and J, are the indices of the represen- 
tation of SU{2). However, deviation from D-flat direction is strongly restricted. Therefore we restrict dynamics of 
scalar fields into the space of flat directions. Hereafter, we consider scalar fields parameterized as 



Hu — 



(16) 



and, 



(17) 



We can easily see that Dy = Di = D2 = = Q and thus Vd by this parameterization and there remains no 
additional physical degrees of freedom. Here F-term in Ref.pJTj vanishes because of D-flat conditions. 

Moreover, we assume 



Alii ~ A// ^ Al22,Al33, 
in other words, neutrino mass matrix has a hierarchy. In this case, potential minima is located at 



\hd\ 
\i>2\ 



(HM/Xlh 
(HM/\hY'^ 



1/2 



l^3l 



0, 



(18) 

(19) 
(20) 
(21) 
(22) 



because 1^2 and acquire heavy masses from hu- Therefore only one slepton and two Higgs have large vevs and other 
scalar fields set in the origin approximately. Hereafter, we adopt a notation, A^n — A^, vi — v and A representing 
the mean value of Al and \h- 



B. CP-violation 



We now comment on CP-violation. In the case only a single flat direction has a nonvanishing vev, there are two 
phase-dependent terms during inflationary era[Tol [T^. 



Am3/2X(p" 



h.c. = 



h.c. 



2|a|i7A|<^|" 
2|A|m3/2A|^|" 



cos{n6tj, + 9a), 



(23) 
(24) 



where <p is the AD field parameterizing the flat direction, 9a and 9a are phases of the coefficients a and A, respectively. 
In this model, (j23p determines the initial condition just after inflation and (|24p determines the eventual potential 
minimum in the angular direction for the following dynamics. Therefore 9a — 9a represents the CP-violation. When 
H ^ '713/2, in other words when ([24p is negligible, there is no CP violation, and lepton asymmetry cannot be induced 
because of Sakharov's condition As decreases to make the ratio of A-term from hidden sector to the total scalar 
potential for AD field larger, CP-violation becomes effective and lepton asymmetry can be generated. In our model, 
on the other hand, we have five phase-dependent terms during inflationary era. 



Xi 



■m3/2ALi''^hl + h. 



2M 
Afl_ 
2M 



HaHHjHl + \i.c. 



^my,AHHjH^ + h.c. 



Al.2> , ,2, 



^WlM^KI^ cos(2(?, + 29u + 9i) 

ALm3/2 

M ' 
M 

AffTO3/2 



Al\\vV\KVcos{29, + 29u + 9l), 
aH\\hu?\hd\''cos{2eu + 29d + eu), 



M 



\AH\\K\''\hd\''cos{29u + 29d + 9H) 



M2 



i>r\hu\'\hd\'cosi29,~29d), 



(25) 
(26) 
(27) 
(28) 
(29) 
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where 9^, Ou, ^d, Qi, 9h, and 0h are phases of h^, hd, a^, A^, an, and Ah, respectively. Four of them are A-terms 
and the last one is from F-term. In inflationary era, (PSj) . and (|29p are dominant and they determine the 
initial condition of scalar fields in the angular direction for the following dynamics. Moreover, in this case, the time 
dependence of (|25p and (|27p change after inflation, on which we comment in the next section, therefore the potential 
minimum for angular direction is varied. In this sense, CP is violated strongly just after the end of inflation. In our 
model, therefore, lepton asymmetry can be induced just after the end of inflation. 



C. Stability of multiple flat directions 



In the remainder of this section we demonstrate the necessity of considering multiple flat directions. When a flat 
direction acquires a large vev, many other flat directions acquire masses from F-term. Consequently such flat directions 
cannot take large values, so that single flat-direction description can be justified. However, there are directions that 
do not acquire F-term from other flat directions within the MSSM [l^. The case of LHu and HuHd is its typical 
example. 

The mass terms of scalar fields are 



neglecting soft terms from hidden sector. Setting the D-flat condition (fTT]) . it reads 

- (c, -f c„)i/>|2 - (c„ + C<i)H>d|^ 



(30) 



(31) 



li Cn + Cu > and c„ -I- Cd < 0, /id is stable at the origin and we conclude that hd falls down to the origin and we 
can take the one-flat direction description. However, if Cj/ -I- c„ > and Cu + Cd > 0, both v and hd are unstable and 
cannot stay at the origin. We can no longer take the description. Moreover, even if only LH^^ direction is selected in 
the beginning, Hd is unstable at the origin for large parameter region. We can see this by evaluating mass matrix for 
Hd around the local minimum of LH^ fiat direction parameterized as (j) , 



H' 



Cd 
Cd 



M 



Re(aH) 



101- 
\hH 



2M 
XhH_ 
2M 



mJ I'^I (,0 1 

(Im(0))2 -2Re(</>)Im(0) 
-2Re(0)Im(0) -(Re(0))2 -l- (Im(0))2 

(Re(0))2 - (Im(0))2 -2Re((?!))Im(0) 
-2Re((?!))Im(0) -(Re(0))2 (Im(<?!)))2 

-2Re(0)Im((/)) -(Re(0))2 + (Im((/)))2 
-(Re((?!)))2 + (Im(0))2 2Re(</>)Im(0) 



The eigen values of this matrix are 



^1 W 



CdH^ ± 



|0|4 + Re(aH; 



2M 



I0P 



Im(aH) 



\hH_ 
2M 



1/2 



(32) 



(33) 



As a consequence, hd is unstable around LH-^ fiat direction with a wide range of parameters. For example, even 
when Cd > 0, sufficiently large an or can make hd unstable. (Note, however, that too large Xh makes i> stable.) 
Therefore, the situation where multiple flat directions have nonvanishing vevs is more natural than the case with 
single flat direction. 



III. DYNAMICS OF SCALAR FIELDS AND LEPTOGENESIS 



In this section, we describe the motion of scalar fields and estimate baryon-to-entropy ratio. Here we consider only 
homogeneous mode and neglect fiuctuations around it because the curvature of the potential is much larger than the 
Hubble parameter and quantum fiuctuation is sufficiently suppressed in this case. 

In the Friedman universe with the spacetime metric, ds"^ = dt^ — a?{t)dx'^ , scalar fields (pa obey the equation of 
motion. 



(34) 
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Parameterizing as ((T6)) . there remains a D-term, 

VD = {g^ + g'^){\v\''-\K\'' + M''f. (35) 

Numerical analysis shows that this potential energy density remains small and that deviation from flat directions 
by this parameterization does not have great influence on the dynamics of scalar fields and the resultant lepton 
asymmetry. This result justifies our definition of the field variables (|16|1 with which all the D-terms except (|35p vanish 
from the beginning. If ^ d'^V{(t)) / d(j)ad(j>*, the friction term is dominant and the fields are overdamped. On the 
other hand, if <C d'^V{(l))/d(f>d(f>* , the inertial term is dominant and the fields are underdamped. 

In inflationary cosmology, the Hubble parameter evolves as follows. In inflationary epoch, the Hubble parameter is 
approximately constant and we call it Hinf. As inflation finishes, the inflaton starts oscillating around its minimum 
and the Hubble parameter decreases as in the matter dominant era, H — 2/3i. Meanwhile the inflaton decays into 
other light particles gradually, and finally the universe turns to radiation dominant. Note that the cosmic temperature 
reaches the highest just after the end of infiation as a result of partial decay of the inflaton even if preheating is not 
operative. As a consequence, we cannot neglect the effect of the high-temperature cosmic plasma on the scalar fields 
even during the inflaton oscillation dominant regime p^. [isl. KE^ . 



A. Inflationary era 

During inflationary epoch, scalar fields go to a minimum of the potential gradually. In our model, the potential ^ 
is very complex and it is nontrivial how many minima our model has and which minimum is selected by those fields. 
If the number of e-fold is large enough, scalar fields fall into one of the potential minima because d'^V{(j)) /d(f>d(f)* ^ 
HinfM/X ^ ^inf the minimum. Moreover, using numerical analysis we have confirmed that the minimum is 
unique except for gauge freedom and phase inversion. Figures [T] show typical example of the value where scalar fields 
fall from various initial values. Setting sneutrino real by U{1) gauge transformation, /i„ and hd are found to be on 
only two minima whose amplitudes are equal and whose phase differences are n. What is important is that the values 
of the fields are fixed both in the angular direction and in the radial direction. In particular, CP-violation just after 
infiation arises from difference between the initial phases of scalar fields and the time-dependent potential minima in 
the angular direction at the inflaton oscillation dominant era. Note that if 0i = 0h + tt, the potential minima during 
and after inflation are equivalent until the beginning of the oscillation of scalar fields around the origin. In this case, 
CP is not violated at large H. The change of potential minima is discussed precisely in the following section. 

As a result, the vevs of scalar fields are fixed on the minima at the end of infiation, and following dynamics of scalar 
fields is independent of initial value at the beginning of inflation because quantum fluctuations are suppressed due to 
the large effective mass at the potential minima. In other words, there remains no pre-inflationary information. This 
feature is the same as in the one-field case. 



B. Inflaton oscillation dominant era 



After infiation the universe enters the inflaton oscillation dominant era. The radial components of motion of scalar 
fields in this epoch are almost the same as in the case of single flat direction. 
In this epoch the equation of motion is 

Ct^a + + ^ = 0. (36) 

Neglecting phase dependent terms, and reparameterizing as 

t^z = ln(tA/), tf^'^H-„), (37) 

<i>a^Xa^C^a/{HM/\Y'\ (38) 

becomes [13 



This equation has a fixed point (fa > 

0, ^ (10/9)1/40(0), (40) 
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FIG. 1: The value of scalar fields at the end of inflation, the real part (horizontal axis) and the imaginary part (vertical axis) 
before (top) and after (bottom) gauge transformation. Red crosses represent D, blue x's represent hu, and green stars represent 
hd. The following parameters are used, c^, = 0.8, c„ = 1.0, Cd = 1.2, M/\l — 1.2 x W^Hi„f, M/\h ~ 1.5 x 10^-ffi„/, 
ql = e™''^, and an = i- Initial values are chosen randomly in the range O.OlHinf ~ 100.0Hi„/. 

where ^i"^ ~ (HM/X)^^^ are minima of the potential. This fixed point presents a solution which traces points that are 
a little bigger than the time-dependent minimum of the potential. The initial values are the minima of the potential 
and different with these fixed points. Therefore scalar fields oscillate around the time-dependent fixed point with 
the period Az ~ 1. Due to the time dependence of the potential minima, the oscillations do not redshift away even 
in the presence of cosmic expansion. As a result, scalar fields decrease with time as (f>a ^ {HM/Xy^^ ~ {M/Xty^^ 
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oscillating around the fixed point. 

When a soft term from hidden sector, \4>a\'^ , or, thermal terms [T^ [Tsj. 



h\<f>a\<T 



or OgOfs (T)^T^ In 



rp2 



(41) 



dominate the Hubble induced mass terms, the minima of the potential are at the origin and scalar fields begin 
oscillating around it. Here k is an index of fields which couples to flat directions, fk is its coupling constant, Ck and 
ttg is numerical factor of order unity, and as{T) — 5s(T)^/47r. {gs is gauge coupling of SU{3). ) 



When scalar fields oscillate around the origin, the amplitudes of scalar fields change with Hubble parameter as 



H (when driven by soft mass term from hidden sector) 
|(^a| oc <( i/Vs (when driven by thermal mass term) 
(when driven by thermal log term). 



(42) 



Figure [2] shows typical time variation of scalar fields. The center of oscillation deviates from (|40|) because the 
complicated mixings of scalar fields change the fixed points Xa- 
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FIG. 2: Time variation of amplitudes of scalar fields. Horizontal axis is logarithmic time and vertical axis is amplitude of scalar 
fields. Red solid line is v, green dashed line is hu, blue dotted line is hd. The following parameters are used, = = = 1.0, 
Hir,f = 1.0 X lO^^GeV, M/\l = 5.0 x lO^^GeV, M/\h = 7.5 x 10^'^GeV, = e'''''^ an = i, and reheating temperature Tr 
is 5.0 X lO^GeV. 

However, the angular motions of scalar fields are different from one-field case. For this purpose, we first comment on 
the Hubble induced A-terms in this era. The Hubble induced terms in inflationary era, —CaH'^\(j)a\'^ and {aXH/M)(f)^ + 
h.c, are from the couplings with \Fj\'^ and Fj, respectively. These couplings exist even in the inflaton oscillation 
dominant era. In this regime, inflaton oscillates around the potential minimum with sufficiently shorter time scale 
than the motions of AD fields. The induced potential, therefore, can be acquired by replacing and Fj with 

their time average. Generally (((5/)^) is dominant in \Fj\'^ because the linear term (SI) vanishes by time averaging, 
where SI is the fluctuation of inflaton from the potential minimum and () represents time average. Moreover, when 
the F-term of inflaton includes second order monomials of the inflaton, is also dominant in As a 

consequence, Hubble induced A-terms have the same dependence on time as the Hubble induced mass terms. Since 
(^/) = (l^/P) — 3H^M^/2 by the Virial theorem, the Hubble induced A-terms can be written in the form. 



(43) 
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where k is a numerical factor that is determined by the inflation model and here we assume that its magnitude ia 
order of unity. For example, if we take the superpotential of the form as a inflation modelp2i]. 

W{I,X) = v^X{l-gI^), (44) 

where / is inflaton, and X is auxiliary field, v is vacuum energy, g is coupling constant whose dimension is —2, the 
important F-term is 

Fx=v^il~gl^). (45) 

The resultant potential is 

y(/, X) = v^[{l - giy + Ag^I^X^] (46) 

In inflationary epoch, the vacuum energy dominates the universe, V = V{0,X) = v"^. Therefore, the Friedman 
equation reads 

v^ = 3HlfM^ (47) 
After inflation, the inflaton oscillates around the potential minimum. As the potential minimum is 

lo = VT/g, (48) 



we expand 



Therefore, the F-term can be written as, 



and we take the time average of them. 



I=y/l/]j + 6I (49) 

Fx = -v^[2^SI + giSI)% (50) 

\Fx\^^^gv'\6I\', (51) 

\{Fx)\^v'g{\5I\'), (52) 

{\Fx\')^4gv\\Sm. (53) 



As a consequence, these two time averages are related as, 



\iFx)\ = {\Fx\')/{Av') = ^^^, (54) 



and we can write the Hubble induced A-term as, 

{Fx)K'''^D^W ~ 0^ + h.c. (55) 

As the amplitude of AD field decreases in proportion to H^^^, the Hubble induced A-terms fall off more rapidly than 
the mass terms. Note that the inflaton X has i?-charge 2 and the Kahler potential of ([T]) type break i?-symmetry. 
Here we assume that i?-symmetry of the inflaton is broken by some mechanism of high energy physics. 

In the one field case, the change of A-term discussed above does not make the change of the minima potential 
in the angular direction at large H and it is fixed until just before the beginning of oscillation around the origin 
because the effect of CP violation is very small at that time. As we commented in the previous section, in our model, 
the phase difference of coupling constant of two A-term, ((23)) and ((24)) yields CP violation. In contrast to the one 
field case, there exist phase dependent terms in the potential ((^^ and the Hubble induced A-terms. Moreover the 
change of A-terms ((55)) induces the change of potential minima in the angular direction. Therefore scalar fields also 
start oscillating around the minimum in the angular direction as in the radial direction at the onset of the inflaton 
oscillation. The equation of motion for the phase, 9a, of a scalar field 4'a, 
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leads 



4 dV 

+ -7-^— = 0. 



9z2 9 

In this case, the relevant part of the potential is of the form of 



V 



sin(2^?, - 29,1). 



(57) 



(58) 



Here we approximate that all the scalar fields have the same amplitude. Therefore, the phases of v and hd oscillate 
around its minimum without damping. The period of its motion is also Az = However (j58p does not depend 

on the phase of hu and the only phase dependent terms, or the Hubble induced A-terms, for hu decreases rapidly 
and the potential for /i„ in the angular direction becomes flat. Therefore angular motion of is fixed immediately. 
When the oscillation around the origin begins, the angular momentum is fixed to its value at that time because phase 
dependent term in the potential becomes negligible and it becomes the constant of motion. Figures [3] show typical 
motion of scalar fields. Scalar fields oscillate not only in the radial direction but also in the angular direction even 
before the oscillation around the origin begins. The angular motions of v and hd are swinging, on the other hand, 
that of hu is not. 




'0.6 '0.4 -0.2 
Real part of V ((HM/X)'"^) 




'1.2 -1 -0.8 -0.6 -0.4 -0.2 
Real part of fi„ ((HM/)l)""^) 




-0.2 -0.1 0.1 0.2 
Real part of fi^j ((HM/i)""^) 



FIG. 3: Motion of scalar fields {v,hu, and hd from left to right) during inflaton oscillation dominant era. The value of fields 
is rescaled as x = (f)/{HM /Xf'^'^ . The following value of parameters are used, = Cu = = 1.0, Hinf = 1.0 x lO^^GeV, 
M/Al = 1.0 X 10^°GeV, M/\h = 6.7 x lO^^GeV, hl = e'"^^, as = i, ^ = 8 and Tr = 5.0 x lO'' GeV. 



C. Lepton asymmetry 

Now we turn to the lepton asymmetry. Lepton- number density is expressed by the time component of the 
Noether current, or the operator, 

nL(t) = i{i>*i)~i)*i>)^2\i>\'^0^. (59) 
Therefore, using equation of motion we can derive the equation, 

^^V \ 

fiL + -iHuL = 2Im —V . (60) 



In the integral representation, wc find 



dv 



^dt) = —r. / dt'a{t'flu, — J> . (61) 



■ / dt'aU'fl^ ( 



a{tf 

In the case of single flat direction, LH^ [l3|, we evaluate the lepton asymmetry as 



"LW = ^y^ dt'a{t'f^^\^\UHAe^ + ea), (62) 
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where (p = D = hu represents AD field. Here we consider only the hidden sector induced A-term and neglect the 
Hubble induced A-term because AD field is located at its minimum in the angular direction and its contribution is 
negligible. The net lepton-to-entropy ratio is [l3| 



TlL TrM ( TO3/2 



s 12AlA/2 \Ho 



^eff. (63) 



Here Seff — sin(4^?0 + 9a) represents an effective CP violating phase. It has a factor 111^/2/ Hose because the gradient 
of the potential to the angular direction is smaller than that to the radial direction by a factor of m^^2/ Hose, which 
represents the dilution of CP violation. In the case of multiple flat directions, on the other hand, 



nL{t) =^ / dt'a{t') 



t 

i\3 



^ -H^\AlM^\K\^ sin(20, + 29^ + 9l) 



nMHinf 

+ ^my2\ALM^K\^ sm{29, + 29^ + 9l) + -^^\i)\'\K\''\hd\'' s\x,{9, - 9^) 

dt'a{t'f^H\aL\\(l>fsm{59)[l + 0{K-') + 0{m3/2/H,nf)]. (64) 

Here we set tf as the time when inflation flnished We also use the fact that |j>| ~ \hu\ — \hd\ and |(/)| represents the 
mean value of them. 

If sneutrino were in its phase minimum, liii[{dV/ di>)i>] ~ 0, lepton asymmetry would not develop. In our model, 
however, sneutrino oscillates around the phase minimum during inflaton oscillation dominant era and does not stay 
at the minimum. The total lepton asymmetry can, then, be estimated as 

4 /■* MH^ 
n,{t)^—^ dt'a{t'f\aL\^—siTi{5e) 

<^[tr Jtf A 

3 A/T pr2 



3U(ty) \''L\^^sm[5,ffitosc))[l + OiHosc/H,nf)], (65) 



where tosc is the time when the oscillation around the origin starts and S9 = 29^ — 29d- Here 5eff{t) is the total 
phase of the lepton asymmetry at t that is due to the oscillation of 59. For t > tosc, since the amplitudes of scalar 
fields decrease more rapidly than H, the contribution to the net lepton asymmetry in this epoch is suppressed by the 
factor 1713/2/ Hinf- Here we neglect the Hubble induced A-term and the soft A-term because their contributions are 
very small. Moreover, we do not have a factor \n{tosc/tf) because of the time-dependence of 69. 

When reheating is finished, lepton-to-entropy ratio n^/s is fixed because both lepton number density and entropy 
decrease in proportion to a^^ afterwards. The final lepton-to-entropy ratio is therefore given by 

s s 

= TTT-TTTx ^Mr^ r^^sm[Seff{tosc)) 

{2TT/Ab)g^sTff \a{tR) J A V V 

- (2./45),..r^ ^i4:;MrJ l«.|^sm(..,,(W)) 
2\aL\ ( .9, A TrM . / N 



where g^, is the total number of effectively massless degrees of freedom, Ir is the time when reheating is finished and 
Mq is reduced Planck mass. Without numerical factor of 0(1), we find. 



^(t)^0.1x|^sin(5e//(W)). (67) 

G 



As a result, the sign of lepton-entropy ratio is determined by the information ait — tosc but its order is determined only 
by the cut-off scale M, the smallest coupling A, and the reheating temperature Tr. Contrary to one-field description, 
there is no dilution factor 171^/2/ Hose, because in the present case CP is violated by the Hubble induced terms and 
is not weakened by thermal effect, which exists in the case of single flat direction. This implies that the net lepton 
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asymmetry is determined almost by the parameters associated with the couphng between inflaton and fiat direction 
and the low-energy SUSY breaking mechanism affects only a rather arbitrary phase factor (5e//(iosc)- 

While the number density of hd is created in the same mechanism, that of is rather small because it is created 
from the vanishing Hubble induced A-tcrms. This is the different feature with Rcf. [17]. 

Figure |4] shows typical time variation of the lepton asymmetry and the number density of Higgs fields. Lepton 
asymmetry and number density of hd par comoving volume fluctuates just after inflation finished, and fixed when 
oscillation around the origin begins. On the other hand, the number density of per comoving volume is fixed 
immediately after infiation because of the vanishing Hubble induced A-terms as discussed above. Figure [5] shows the 
net dependence of lepton asymmetry to phase difference between and an- When exactly — an ~ 0, or tt, the 
net lepton asymmetry is equal to zero, as it should be. 
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FIG. 4: Time evolution of lepton asymmetry and number density of Higgs fields par comoving volume. Horizontal axis is 
logarithmic time and vertical axis is lepton number density par comoving volume divided by entropy par comoving volume at 
the time when reheating is finished. Red solid line is lepton asymmetry, green dashed line is number density of hu and blue 
dotted line is that of hd- The following parameters are used. c„ = Cu ~ Cd — 1.0, Hi„f = 1.0 x lO^^GeV, M/Xl = 1.0 x lO^^GeV, 
M/Xh = 1.5 X lO^^GeV, ul = e'^"''^ an = i, and Tr = 5.0 x lO^GeV. 

The sphaleron effect converts lepton asymmetry to baryon asymmetry as [2^ 

^ = 4^. (68) 
s 23 s 

Therefore present baryon-to-entropy ratio is represented as a function of reheating temperature and lightest neutrino 
mass, 

— ~ -0.02 X ^— ^sm(5e//(W)) 
Here we have taken the phase factor sin(5e//(^osc)) = 0(0.1) and used the equation ^ assuming tan/3 — 0(1). 



IV. THE CASE OF NO HUBBLE INDUCED A-TERMS 



In the discussions so far, we have concentrated on the case with the Hubble induced A-terms. As mentioned in 
section fll Al however, there is a case where Hubble induced A-terms is forbidden by some symmetry. In this section, 
we discuss the infiuence of no Hubble induced A-terms and comment on its cosmological implication. 
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phase difference between and a|_| (rad) 



FIG. 5: The dependence of lepton asymmetry to phase difference between and an- Horizontal axis is phase difference and 
vertical axis is the net lepton-to-entropy ratio. The following parameters are used. c„ — Cu — Cd — 1.0, Hi„f — 1.0 x lO^^GeV, 
M/Xl = 1.0 X 10^°GeV, M/\„ = 1.5 x lO^^GeV, and Tr = lO^GeV. 



A. CP-violation and stability of a flat direction 

In section |lTl we have seen the CP-violation and stabiUty of flat direction in our model. If the Hubble induced 
A-ternis ([25]) and (P7|) are forbidden, the situation is changed. During inflation only is the phase dependent term. 
The initial conditions of scalar fields in the angular direction for the following dynamics are not determined completely. 
The hidden sector induced A-terms ([26]) and ([28]) are responsible for CP-violation. The degree of CP-violation has 
indcfinitcness because of arbitrariness of initial phase. Furthermore, lepton asymmetry is induced at lower energy. 

As for the stability of a flat direction, we can set ai, — an = in the mass matrix (|32)) and the eigenvalues (13311 . In 
(155)1 the large coupling constants of the Hubble induced A-terms are important for the instability of a flat direction. 
Therefore, the parameter space where we must consider multiple flat directions is smaller than in the case of existence 
of Hubble induced A-terms. Hereafter, however, we discuss the case where we must consider multiple flat directions. 



B. Dynamics of scalar fields 

We turn to the difference of dynamics of scalar fields. While the dynamics in the radial direction is not changed, 
the dynamics in the angular direction is different. 



1. Inflationary era 



During the inflationary era, the relevant potential for the angular direction is the cross coupling in the F-term, 

^^H'K\^\hd\^cosi29, - 29d) e V. (70) 

We can easily see the potential minimum for the phase of v and is 

0^-0d = {2n+l)Tr (neZ) (71) 

However, the potential for the phase of is flat and the initial phase of for the following dynamics is not determined 
by the potential. (Strictly speaking, considering the gauge freedom, the true physical degree of freedom is the phase 
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difference between /i„ and 9u — 9^.) Therefore, the initial values for the following dynamics have indefiniteness. 
In other words, there remain pre-inflationary informations. In the following section, we see the influence of the 
indefiniteness for the resultant lepton asymmetry. 

2. Inflaton oscillation dominant era 

We have seen that scalar fields begin to move in the angular direction just after the inflation when the Hubble 
induced A-terms exist. The source of this feature is the difference between the potential minima in the angular 
direction during inflation and at the inflaton oscillation dominant era. When the Hubble induced A-terms do not 
exist, however, the potential minima in the angular direction are not different between during and after inflation 
because (fTOl) is the only phase dependent potential for large H. In this case, therefore, scalar fields begin to move 
in the angular direction when the oscillations around the origin start as is the case of one field [l^]- We can see this 
feature in Fig. [6l 




FIG. 6: Motion of scalar fields during inflaton oscillation dominant era in the complex plane, D, hu, hd from left to right. The 
value of fields is rescaled as x = 4>/ {H M / XY^'^ . The following parameters are used. Cv = Cu = Cd = 1.0, Hinj — 1.0 x lO^^GeV, 
M/Xl = 1.0 X lO^OQeV, M/Xh = 3.3 x 10'^°GeV, Al = e'"''^ Ah = i, and Tr = lO^'GeV. 



3. Lepton asymmetry 

Now we comment on the resultant lepton asymmetry. We can estimate it by eq. (|64[) . From this equation, we can 
see that the deviation from the potential minima in the angular direction is essential for generating lepton asymmetry. 
While deviation from potential minima exists just after the end of inflation when the Hubble induced A-terms exist, 
it does not exist until the beginning of the oscillation of scalar fields around the ori gin when the terms do not exist. 
In this case, moreover, the integrand is practically the same as that of one field case [lOl [l5|. Therefore, the resultant 
lepton asymmetry is [l5|, 

UL MTr fmy^\ 

where 5eff is the effective CP- violating phase. We can see that the suppression factor 777.3/2 /-ffosc is restored. Figure 
[7] shows the typical time variation of the lepton asymmetry and the number density of Higgs fields. This feature is 
the same as one field case [lo|, [Hi • 

The CP-violating phase 5eff depends not only on the couplings of A-terms, and Ah, but also on the initial 
phase of h^, which represents spontaneous CP-violation, as mentioned above. Figure [5] shows the dependence on the 
initial phase of hu- 

C. Baryon isocurvature perturbation 

Now we comment on baryon isocurvature perturbation 25). During inflation, in this case, the phase of hu, 0^, is 
effectively a massless field while other physical degrees of freedom have masses of order of H . Therefore, we must 
consider baryon isocurvature perturbation from the quantum fluctuation of 9u as is the case with Ref. pjj . 
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FIG. 7: Time evolution of lepton asymmetry and number density of Higgs fields par comoving volume. Horizontal axis is 
logarithmic time and vertical axis is lepton number density par comoving volume divided by entropy par comoving volume at 
the time when reheating is finished. Red solid line is lepton asymmetry, green dashed line is number density of and blue 
dotted line is that of hd- The following parameters are used. c,y = Cu ~ Cd — 1.0, Hi„f — 1.0 x lO^^GeV, M/Xl = 1.0 x lO^^GeV, 
M/X„ = 3.3 X lO^^GeV, Al = 1.0, Ah = i, and Tr = 1.0 x lO^GeV. 




initial pliase of 11^ [rad] 



FIG. 8: The dependence of lepton asymmetry to the initial phase of Horizontal axis is the initial phase and vertical axis is 
net lepton entropy ratio. The following parameters are used. c„ = Cu = Cd = 1.0, Htnf = 1.0 x lO^^GeV, M/Xl ~ 1.0 x lO^^GeV, 
M/Xh = 3.3 x lO^^GeV, and Tr = lO^GeV. 
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Assuming other degrees of freedom set at their potential minima, the equation of motion of 9^ is 

On + iHOu - ^A0„ = 0, (73) 

or we set (p„ = \hu\9u-, 

(Pu + iHifu - ^Av5u = 0, (74) 

Here we have taken the spatial derivative into account. Note that this equation of motion is applicable only if the 
fluctuation of 9u is small enough. If we decompose the quantum fluctuation of ipu, 0u, as 



/ d rZ 

= / A^ZvU, + 4<fcWe""1, (75) 



(27)3/2 

^uk(t) satisfies the equation 

ifiukit) + iHifiukit) + k\-^ipuk{t) = (76) 
Here Ofc and at are annihilation and creation operators, respectively. In the de Sitter space, the solution of (|76[) is 



^uk{t) - ^H,nf{-v)'^''H^yli-kv), V = -iH.nfait))-\ (77) 

where Hy]^ is the Hankel function of the first kind. Here the vacuum state is chosen so as to coincide with the 
usual Minkowski vacuum at ry — > — oo. Application of ([75|l or (|74p is therefore justified since the fluctuation of 9u, 

9u ~ H,r,f/\K\ ~ [H,r,f\/Mf/^ « 1. 

Using the following expansion of the Hankel function, 

i/«(-M^-^(^) for-fc,y«l (78) 

we can see the power spectrum of the perturbation of 'P{k) = [kr' /2'n:'^){if*^^Lpuk) , is 

(79) 

The power spectrum of the fluctuation of the phase of /i„ is, therefore. 



3 



2 



Vs(k).^J.W^y^^^ .y^y (80) 

The fluctuation of the phase of /i„ generates the fluctuation of the net lepton asymmetry, SriL^k), which we can see 
from Fig. [51 

+ SuLik) ^ - i^j \aL\^^siniSeffitosc) + 2Uk))- (81) 



Therefore baryon isocurvature perturbation is 

StiLik) 



We can see from the observational constraint ub/s ~ 0.9 x 10""'^'' [l| and 



2 COt{de ff{tosc))9uik). (82) 



c; \ A , r2 

sin(5,//(W))=^-^;^xlO-^ (83) 
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The power spectrum of baryon isocurvature perturbation is 



1 



I 7:VM ] I A2M4 j 



(84) 



Therefore larger reheating temperature and cut-off scale make larger isocurvature perturbation. Note that Eq. 
means that 5AM| x 10^ < TrM. 

The constraint on the baryon isocurvature perturbation in terms of the ratio between the power spectrum of 
matter isocurvature perturbation Vsik) and that of curvature perturbation V-jzik) is known as [2y| Ba < 0.31. Here 
Ba = \/'Ps{k)/'Pn{k). Using the observational data, Vn = 2.4 x 10"^ [l|, and the relation, 

(^^^VbW, (85) 
we find that parameters of this mechanism are constrained as 



cot^ (5,ff{t,,,)) - X [-^^-y^ 1 j < 7.0 X 10-«. (86) 



M 



Here VLjy and U,,n are the density parameters of baryon and matter, respectively, with their observed ratio f2b/r2,„ = 
This constraint requires small Hubble parameter during inflation. Moreover, small cot ^(5e// (iosc )^ , or small 
reheating temperature or small cut-off scale is required. 



V. CONCLUSION 



In this paper, we have obtained three important results about Affleck-Dine leptogenesis via multiple flat directions 
with non-renormalizablc superpotential and vanishing renormalizable F-terms. First, when multiple flat directions 
have negative Hubble induced masses, we can no longer parameterize flat directions in terms of one complex scalar 
field and multi-dimensional motion of scalar fields must be considered. The scalar potential, however, has unique 
minimum except for gauge freedom and phase inversion if the Hubble induced A-terms exist. Therefore the degree of 
the mixing of flat directions is determined only by the shape of the potential without ambiguities, and the scalar-field 
dynamics in the post inflationary era are deterministic. Thus in the AD mechanism via multiple flat directions with 
nonrenormalizable potential to realize the nonvanishing field values, there remains no pre-infiationary information 
provided that infiation lasts long enough. In the case where the Hubble induced A-terms are absent, on the other 
hand, one of the physical phases of scalar fields takes an indefinitive value and its relative phase to the eventual 
potential valley acts as a source of CP-violation as in the case of single-field without the Hubble induced A-term. 

Second, there are CP-violation terms even for large H due to cross coupling of scalar fields in non-renormalizable 
F-terms and the Hubble induced A-terms. Although the Hubble induced A-terms decreases rapidly after the end of 
infiation, they can give angular momentum to scalar fields. Therefore lepton asymmetry is generated just after the 
end of inflation. In particular, there is no suppression due to thermal effect [l3|- Net lepton-to-entropy ratio does not 
have suppression factor m^/2l Hose, and the low-energy SUSY breaking mechanism, in particular, the gravitino mass, 
has practically no influence to the net lepton asymmetry. Since the Hubble induced A-terms play the critical role, 
the absence of them destroys this feature and the suppression factor m-^iij Hose is restored. 

Finally, in the case without the Hubble induced A-terms, this mechanism generates the baryon isocurvature pertur- 
bation and observational constraints narrow the allowed parameter space. In particular too large Hubble parameter 
during inflation is disfavored by this constraint [25| . 
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